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Newtonian  fluid  and  as  a visco-plastic  material. 
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defined  and  an  order  of  magnitude  analysis  is  performed  for  both  the  Newtonian 
fluid  and  visco-plastic  representation  of  the  stress  tensor.  The  two  models  for 
the  stress  tensor  are  compared  and  shown  to  be  equivalent  for  small  values  of  the 
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I . IkTTRODUCTION 


A set  of  partial  differential  equations  is  formulated  to  describe 
the  shaped  charge  jet  behavior.  Somewhat  similar  equations  are  given  in 
the  various  hydrocodes  but  these  codes  employ  various  assumptions  re- 
garding the  modeling  of  the  stress  tensor  and  use  mixed  Eulerian- 
Lagrangian  systems.  However,  hydrocode  analyses  include  compressibility 
effects  using  semi -empirical  equations  of  state  relating  the  density, 
pressure,  and  internal  energy. 

In  this  report,  incompressible  flow  is  assumed  at  the  onset,  since 
one-dimensional  hydrodynamic  analyses  indicate  that  compressibility 
effects  are  not  dominant  during  the  jet  development  or  jet  penetration 
process  (see  for  exan^le  Allison-VitaliS  DiPersio-Simon-Merendino^, 
Eichelberger^,  and  Chou-Carleone-Tan2io-Cicarelli‘*) . The  major  influence 
of  compressibility  is  in  the  liner  collapse  process  as  modeled  by  Chou- 
Car leone- Karpp^®  especially  in  regard  to  a jet  - no  jet  criterion. 
However,  once  it  has  been  established  that  a given  liner  geometry  and 
liner  material  will  form  a coherent  jet,  the  collapse  process  can  be 
satisfactorily  modeled  as  an  incompressible  flow  as  done  by  Pugh- 
Eichelberger-Rostoker^ . 

In  this  report,  the  one-dimensional  hydrodynamic  equations  are 
extended  to  an  axisymmetric  hydrodynamic  theory  using  a Newtonian  fluid 
model  for  the  stress  tensor.  Also,  the  material  strength  effects  are 


E.  AViieon  and  R,  Vitatx,  "A  New  Method  of  Computing  Penetration 
Variabtee  for  Shaped-Charge  Jets,"  BRT  Report  No.  1184,  January  1963. 
rAD  #40048t) 

^R.  DiPereio,  J,  Simon,  A.  Merendino,  "Penetration  of  Shaped-Charge  Jete 
Into  Metallic  Targets,"  BRL  Report  No.  1296,  September  1966.  (AD  #476717) 

^Eiohelberger,  R.  J.,  "Re-examination  of  the  Theories  of  Jet  Formation 
and  Target  Penetration  by  Lined  Cavity  Chargee, " Carnegie  Institute  of 
Technology,  CEL  Report  No.  1,  June  1954. 

**?.  C.  Chou,  J.  Carleone,  C.  A.  Tanzio,  and  R.  D.  Cioarelli,  "Shaped 
Charge  Jet  Breakup  Studies  Using  Radiograph  Measurement  and  Surface 
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^P.C.  Chou,  J.  Carleone,  and  R.  Karpp,  "The  Effect  of  Compressibility 
on  the  Formation  of  Sf^ed  Charge  Jets,"  Proceedings  of  the  First 
International  Symposium  on  Ballistics,  Orlando,  Fla.,  Nov.  13-15,  1974, 
sponsored  by  the  ADPA. 

C.  Chou,  J.  Carleone,  and  R.  Karpp,  "Criteria  for  Jet  Formation 
from  Impinging  Shells  and  Plates,"  J.  of  Avplied  Physios,  Vol.  47, 

No.  7,  July  1976. 

"^E.  M.  Pugh,  R.  J,  Eiohelberger,  and  N.  Rostoker,  "Theory  of  Jet 
Formation  by  Charges  with  Lined  Conical  Cavities,"  J.  of  Applied  Physios, 
Vol.  23,  No.  5,  May  1952. 


included  directly  for  the  stress  tensor  modeled  as  a visco-plastic 
material.  In  this  sense,  the  incompressible  Newtonian  fluid  or  visco- 
plastic model  represents  an  extension  of  the  incompressible  one- 
dimensional theories  of  Allison-Vitali^ , or  DiPersio-Simon-Merendino^ . 
Alternately,  this  model  represents  a relaxation  of  the  hydrocodes  to  an 
incompressible,  but  viscous,  flow  representation. 

The  Newtonian  fluid  model  is  assumed  since  for  low  strain  rates,  the 
viscosity  of  some  metals  depends  weakly  on  the  strain  rate,  approxi- 
mately justifying  a Newtonian  fluid  model®.  The  Russians,  notably 
Godunov- Deribas -Mali®  and  Mali-Pai-Skovpin®  consider  the  viscous  proper- 
ties of  the  jet  metal  to  be  important. 


Once  the  governing  equations  have  been  established,  they  will  be 
non-dimensionalized  and  an  order  of  magnitude  analysis  will  be  performed 
to  further  simplify  the  final  set  of  equations.  This  simplified  set  of 
equations  may  admit  approximate  analytical  solutions.  However, 
determination  of  the  appropriate  boundary  conditions  pose  the  major 
difficulty  in  obtaining  an  approximate  solution  or  an  accurate  numerical 
solution. 

The  same  basic  equations  may  be  applied  to  the  jet-target  inter- 
action where  the  system  of  equations  apply  independently  to  both  the 
free  jet  and  the  target  and  are  coupled  by  appropriate  boundary  condi- 
tions at  the  jet- target  interface.  In  this  case,  a function  of  the 
ratio  of  the  jet  and  target  densities  is  important,  but  incompressible 
flow  may  still  be  assumed^ *^.  The  viscous  terms  are  also  important  in 
the  jet-target  interaction. 

The  remainder  of  this  report  deals  with  the  development  of  the  free 
jet.  Specification  of  the  jet  parameters,  in  analytical  form,  prior  to 
impact  with  the  target  will  provide  the  first  step  in  analyzing  the 
target  penetration.  The  jet  formation  will  be  coupled  to  the  liner 
collapse  process  through  initial  conditions  and  boundary  conditions. 

In  turn,  the  free  jet  solution  will  provide  initial  and  boundary  con- 
ditions for  the  target  penetration  problem.  Thus,  the  overall  problem 
of  the  collapse,  development,  and  penetration  of  the  jet  may  be  con- 
sidered as  three  distinct  regions:  liner  collapse;  jet  development  and 
growth;  and  jet  penetration. 


^S.  K.  Godunov i A.  A.  Deribas  and  V.  I.  Malij  "Influence  of  Material 
Viscosity  on  the  Jet  Formation  Process  During  Collision  of  Metal 
Plates Novosibirsk.  Translated  from  Fizika  Goreniya  i Vzryva, 

Vol.  11,  No.  1,  Jan  - Feb  1975. 

®’/.  I.  Mali,  V.  V.  Pat  and  A.  I.  Skovpin,  "Investigation  of  the  Breakdowi 
of  Flat  Jets,"  Novosibirsk.  Translated  from  Fizika  Goreniya  i Vzryva, 
Vol.  10,  No.  5,  Sept  - Oat  1974. 
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The  following  sections  illustrate  that  the  Newtonian  fluid  and 
visco-plastic  stress  tensor  models  are  functionally  identical  for  small 
values  of  the  yield  stress  in  pure  shear  (i.e.,  the  visco-plastic  model 
relaxes  to  a Newtonian  fluid  type  model  for  small  yield  stresses).  Also, 
ordering  analyses  are  applicable  to  shaped  charge  jet  studies  and 
pertinent  non-dimensional  terms  can  be  introduced.  In  addition,  an 
incompressible  formulation  of  the  governing  equations  is  presented. 


II.  BASIC  EQUATIONS  - NEWTONIAN  FLUID 

The  Newtonian  fluid  model  is  considered  first. 

The  basic  equations  governing  the  behavior  of  an  axisynmetric  jet 
are  given  below.  These  equations  are:  the  continuity  equatioii;  Navler- 
Stokes  equations;  and  the  energy  equation.  These  equations  are  formu- 
lated for  an  axi symmetric,  incompressible,  isotropic  Newtonian  fluid 
without  body  forces.  An  Eulerian  coordinate  system  is  employed  with 
the  origin  located  at  the  position  where  the  free  jet  is  assumed  to 
begin  its  formation,  i.e.,  beyond  the  region  of  influence  of  the  liner 
collapse  process. 

The  basic  equations  are: 


1.  Continuity: 


2.  Radial  Momentum: 


i^(rV)  0 , 


^ ^ 3V\  ap  ^ ^ i av 

' ar 


v_  bVi 
,zV 


3.  Axial  Momentum: 


(2  2 
a^u  ^ 1 8u  , a u 

XJ 


4 .  Energy : 

Neglecting  radiation  heat  flux,  allowing  for  internal  heat 
generation,  and  where 


# ■ u 
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the  energy  equation  becomes 


F- 


i 


! 


P 


De 

Dt 


* + 


1 aT  a^T\ 

' ' BzV  ■ 


(4D 


The  nomenclature  is  standard.  Basically  V represents  the  radial 
velocity,  u the  axial  velocity,  and  r,  z are  the  radial  and  axial  co- 
ordinates, respectively.  The  internal  heat  generation  per  unit  volume 
is  denoted  by  Q,  p is  the  dynamic  viscosity  of  the  media  encompassing 
the  jet,  i.e.,  the  air  in  the  case  of  a free  jet,  or  the  target  in  the 
case  of  jet-target  interactions  (i.e.,  penetration^ . Alternately,  this 
viscosity  could  be  interpreted  to  represent  the  viscosity  of  the  metallic 
jet  resulting  from  the  velocity  gradient  within  the  jet.  This  inter- 
pretation follows  from  the  analysis  of  the  viscous  forces  acting  along 
the  free  surface  of  the  metallic  jet®.  Also,  X = - 2.  y by  Stokes 
hypothesis,  p is  the  density,  T is  the  temperature,  '^p  is  the  pressure, 
e is  the  internal  energy  and  k is  the  thermal  conductivity.  The  energy 
equation  contains  temperature  dependent  terms  since  energy  dissipation 
due  to  heat  conduction  was  included  in  the  general  energy  equation  (4) . 


Finally,  two  state  equations  are  required  to  interrelate  the  four 
thermodynamic  variables  e,  P,  p and  T,  e.g., 

e = e(p,  T)  (5) 


and 


P = P(P.  T). 


(6) 


The  system  of  six  equations  involves  six  unknowns  namely,  V,  u,  P, 
p,  e and  T.  For  incompressible  flow,  with  p a known  constant,  one  state 
equation  can  be  eliminated.  Also,  if  the  heat  conduction  terms  are 
ignored  in  the  energy  equation,  one  state  equation  may  be  eliminated. 

In  fact,  for  incompressible  flow  as  assumed  in  Equations  (1)  through 
(4),  the  energy  equation  and  both  state  equations  are  decoupled  from  the 
continuity  and  momentum  equations.  In  this  case,  the  continuity  equa- 
tion and  the  two  momentum  equations  may  be  solved  for  the  three  unknowns 
V,  u,  and  P,  for  a known  p. 


i 
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III.  NORMALIZED  EQUATIONS  - NEWTONIAN  FLUID 

Next,  the  basic  equations  will  be  normalized  and  an  order  of 
magnitude  analysis  will  be  performed.  The  characteristic  dimensions 
will  be: 

L = characteristic  length  = length  of  continuous  jet,  prior 

to  breakup, for  application  to 
shaped  charge  jets. 
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u ^ = characteristic  velocity  = = jet  tip  velocity  for 

application  to  shaped  charge 
jets. 


t„  » characteristic  time  » L/ii 
c max 


Now,  R is  taken  to  be  the  maximum  radial  dimension,  for  example,  the 
jet  radius  (or  the  maximum  hole  radius  when  shaped  charge  jet  penetra- 
tion into  a target  is  considered) . The  equations  derived  by  the  order 
of  magnitude  analysis  are  applicable  to  the  free  jet  region,  exclusive 
of  the  liner  collapse  process. 

It  is  assumed  that: 

-j-«  1,  or  the  maximum  radial  dimension  is  much  less  than  the 
continuous  jet  length; 

1 or  -j-  = r is  assumed  small,  or  for  ordering  purposes, 

R £ r « 1 and  r = OCR),  i.e.,  R is  the  "ordering  symbol"  and  is  assumed 
to  be  small; 

^ = I = 0(1); 

and  define 

— V — u P — 

V = u = — , P = j , and  t * t/t^,  where  u^  is  the  jet 

'^o  ^o  0 u 

" 0 

tip  (or  maximum)  velocity. 


Using  these  characteristic  dimensions, the  continuity,  momentum,  and 
energy  equations  in  nondimensional  form  become: 

i ^ (r  V)  + ^ = 0 , (7) 

r 9r  3z 
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where  Re  “ — = Reynolds  number, 
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\ here. 


e ••=  e/u 


q = Q/p  = = internal  heat  generation  per  unit 

mass, 

C = specific  heat  at  constant  pressure, 


0 * 


T - T. 
1 

AT 


Tj^  = initial  jet  temperature. 


and  AT  * Tq  - T. , where  T is  the  total  temperature  of  the  jet,  and  AT 

u/ 

may  be  taken  to  be  analogous  to  Schlichting^®. 


The  following  nondimensional  groups  can  be  identified: 
U _ 1 


p L u Re 
o 


= inverse  Reynolds  number. 


k AT  1 1 1 

^ „ 3 ■ Re  E P ’ 

p L u^  r 

C AT 

— y = E = Eckert  number,  which  is  related  to  the  frictional 


u 


and  = Prandtl  number  = y . 


heat  and  the  heat  due  to  compression. 


Cp 


H.,  Boundary  Layer  Theory,  6th  Edition,  MaGrau-Hilt, 


1968, 
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At  this  point,  the  conduction  heat  transfer  term  can  be  dropped 
from  the  energy  equation  since  the  shaped  charge  jet  forms  and  develops 
in  a matter  of  a few  hundred  microseconds,  at  most.  Thus,  sufficient 
time  is  not  available  to  conduct  (or  radiate)  heat  in  or  out  of  the  jet. 


The  final  form  of  the  energy  equation  becomes 
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IV.  ORDER  OF  MAGNITUDE  ANALYSIS  - NEWTONIAN  FLUID 


The  nondimensional  equations  are  ordered  as  follows: 


z = 0(1) 

r = 0(R),  where  R is  small, 


0(1) 

0(i) 


Assume 


•(!)-  •© 


or  the  radial  and  axial  velocity  gradients 


are  of  the  same  order  of  magnitude,  as  implied  by  the  co.itinuity 
equation.  Therefore,  the  continuity  equation  is  of  order  one.  The 
radial  moiuentum  equation  becomes 


0(R) 
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where  the  order  of  each  term  is  given  above  it.  Thus,  — =■  is 

negligible,  and  for  the  viscous  terms  to  be  included  the  Reynolds  number 

1 2 

is  assumed  large,  or  = 0(R  ).  Thus, 

I (13, 

at  ar  at  ar  \ar^  r ar  t‘/ 

where  the  pressure  gradient  has  not  been  assigned  an  order  of  magnitude, 
but  each  remaining  term  in  the  equation  is  of  order  R. 

The  axial  momentum  equation  becomes. 


0(1) 

0(1) 

0(1) 

0(1) 

V ^ + 

- 3u 
u — = 

3P  ^ 1 

3^'u^ 

3t 

3r 
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■ 3l  ' 

r 3r 

3Z^  , 

where  can  be  neglected,  and  each  term  in  the  equation  is  then  of 
3z 

order  one,  except  for  the  pressure  gradient  which  has  not  been  considered. 

3P 

From  the  radial  momentum  equation,  — is  assumed  of  order  0(R)  or  the 

2 

pressure  gradient  across  the  jet  is  of  order  R . Thus,  the  pressure 
gradient  across  the  jet  is  small  and  since  the  flow  velocity  outside  the 
jet  is  nearly  zero,  especially  if  the  jet-atmosphere  interaction  is  neg- 

lected.  The  — term  is  also  small, 

3z 

For  the  energy  equation,  further  assume  e = 0(1),  and  q = 0(1),  or 
there  exists  no  a priori  reason  to  suppose  that  the  internal  energy  or 
internal  heat  generation  terms  are  negligible.  Then  the  energy 
equation  becomes: 
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or 


assuming  E » 0(1),  and  neglecting  terms  of  order  less  than  one. 
The  final  equations,  returning  to  dimensional  form,  become 


I J_ 
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The  first  and  third  equations  are  of  order  one.  The  second  equation,  the 
radial  momentum  equation,  is  of  order  R.  The  energy  equation  is  of  order 
one,  but  is  not  required  because  of  the  incompressible  flow  assumption. 
Since  the  pressure  gradient  terms  have  been  previously  argued  to  be 
small,  the  first  and  third  equations  suffice  for  the  determination  of 
V and  u. 

The  final  order  of  magnitude  equations  are  equivalent  to  the 
Prandtl  Boundary  Layer  equations  in  cylindrical  coordinates  and  for 
transient  flow.  These  equations  are  parabolic.  The  original  Navier- 
Stokes  equations  are  elliptic.  An  Eulerian  formulation  has  been 


employed. 

The 

final  equations  account  for: 

1, 

Two-dimensionality 

(axis-symmetry) , 

2. 

Transient  effects. 

and 

3. 

Viscous  effects. 
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However,  an  ideal  Newtonian  incompressible  fluid  is  assuned.  This 
assumption  is  in  contrast  to  that  of  Riney^^,  where  a compressible,  but 
inviscid  flow  was  assumed.  Following  the  one-dimensional  analysis  of 
DiPersio-Simon-Merendino^ , incompressible  flow  is  assumed,  but  the 
second  order  terms  in  the  Navier-Stokes  equations,  i.e.,  the  viscous 
terms  are  retained.  The  viscous  terms  are  retained  for  two  reasons; 
they  represent  dominant  or  high  order  terms  in  the  governing  equations; 
and  these  terms  can  be  modified  (by  removing  the  Newtonian  fluid  assump- 
tion and  Stokes  hypothesis)  to  allow  visco-plastic  stress  models  to  be 
considered.  Also  viscous  effects  are  important  in  the  penetration 
process • 

Note  that  if  steady  state,  inviscid,  one-dimensional  flow  is 
asstimed  the  governing  momentum  equations  reduce  to 


du 

dz 


or 

P 

const , 

2 P 

which  is  Bernoulli's  equation,  which  was  utilized  for  one  dimensional 
jet  studies^ > 2*3. 

Also  note  that  the  so-called  pressure  gradient  term  can  be  modified 
by  assuming  that  the  stress  tensor  can  be  divided  into  three  parts. 
First,  the  static  strength  of  the  jet  is  extracted  from  the  stress 
tensor  and  designated  as  a.  The  remainder  of  the  stress  tensor  term  is 
divided  into  the  pressure  term  and  the  deviatoric  stress  term,  which  is 
then  modeled  as  a Newtonian  fluid.  Thus,  the  axial  momentum  equation 
becomes 


Du 


3P  ^ /afu  1 
' dz  ~ dz  * ^ r 3r/  ' 


or 


Du 
^ Dt 


jj  ^ 3u\ 
,2  r dr)- 


^^Riney,  T,  D.,  "Numerical  Evaluation  of  Hypervelooity  Impaat  Phenomenat" 
appearing  in  High  Velooitu  Impact  Phenomena,  Chapter  V,  Edited  by 
E.  KinaloWi  Academia  Press,  1970. 
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Thus,  the  normal  stress  consists  of  the  hydrostatic  pressure  plus  the 
static  strength  of  the  material.  Again  assuming  a one-dimensional, 
steady  steady,  inviscid  flow. 


L I 

i,  ; 


[■  ■ 

i ■ 

i:  , 


i 

t 

i- 


L 


1 


or 


2 dz  dz 


0 


a + 


P = 


constant . 


This  equation  is  Identical  to  the  Bernoulli  equation  modified  to  in- 
clude strength  effects  as  given  by  Eichelberger^ . Using  this  formulation, 
the  governing  equations  can  be  written  as 


3 (rV)  ^ ^ 
r 3r  3z 


(19) 


^ 1,  3V  , 3V 

+ V t — + U — 

3t  3r  3z 


and 


3t  3r 


3u 
^ ~dz 


e! 

3r 


^ /3^V  ^ 1 3V  V \ 


3P’ 

3z 


+ V 


/3^u  . ]_ 

w "'V 


(20) 


(21) 


where  v = u/p  and  p'  = P + o. 

When  V is  assumed  known,  there  are  three  equations  in  the  three  unknowns 
V,  u and  P'.  The  first  and  third  equations  are  of  dominant  order. 


The  derivation  of  those  equations  assumed  that  the  Reynolds  number 
is  large.  For  hyperyelocity  flow  through  air,  this  is  a valid  assump- 
tion. Appendix  A investigates  this  assumption  for  jet-target  inter- 
action, where  the  target  material  represents  the  penetrating  media  or 
where  the  jet  viscosity  is  used  in  the  free  jet  formulation. 


V.  BASIC  EQUATIONS  - VISCO-PUSTIC  MATERIAL 

For  a visco-plastic  material  the  basic  equations  are: 

13.,,^  3u  f. 

F ^ * 3F  “ 0 ' 


(22) 
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I 


i 

[ 


This  system  has  two  unknowns  in  three  equations.  However,  Equation  (23) 
is  not  a dominant  order  equation,  as  will  be  shown. 

Now  from  Cristescu^^,  for  a visco-plastic  material. 


a 


rr 


0 

zz 


2n* 


du 
3z  ’ 


rz 


zr 


n* 


and 


where  for  the  jet  material, 

n is  the  viscosity  coefficient, 
n*  = variable  viscosity  coefficient,  and 
k * yield  stress  in  pure  shear. 

VI.  NORMALIZED  AND  ORDERED  EQUATIONS  - VISCO-PLASTIC 

MATERIAL 

Now  for  V *■  V/u  , u = u/u  , 

o 0 

r ■ r/L, 

^^CrietesaUt  N.,  Dunamia  Flastiait^i,  Wiley  and  Sons^  N.Y.,  1987^ 
Chapter  X, 
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iniUmtkii 


z 


fl 


z/L, 
t - t u^/L, 

k a 2 » 

pu  L 

Re  a -2-  , 
n 

the  normalized  continuity  and  momentum  equations  become: 
Continuity: 

V)  ^ ^ = 0 : 0(1) 

T 3r  3z 


Axial  Momentum: 


J7 


analysis  and  algebraic  grouping  of  the  terms,  the  right  hand  side  of  the 
above  equation  becomes 

I 

I 

1 


0(1)  0(1)  0(1/R) 


Thus,  the  dominant  terms  of  the  right  hand  side  of  the  axial 

1 '^rz'^ 

momentum  equation  is  governed  by  the  — — ^ terms  and  not  the 


30 

or  normal  stress  terms.  Thus, 

d Z' 
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retaining  terms  of  order  one  and  higher,  or 


P 


+ V 


“ * r 3r 


For  k = 0(R),  all  terms  in  the  final  axial  momentum  equatioii  are  of 
order  one.  If  IT  = O(R^) 


3u  . 

- 3u 

— + 

V — + 

u — 

3t 

3r 

3z 

1 3ul 


_ 2 

or  the  axial  momentum  equation  for  visco-plastic  flow  for  k = 0(R  ), 
i.e.,  small,  reduces  (o  that  of  a Newtonian  fluid,  except  the  viscosity 
term  in  the  Reynolds  iiiiinbcrs  may  differ  between  the  Newtonian  fluid  and 
the  visco-plastic  model  unless  the  viscous  effects  in  the  Newtonian 
fluid  model  are  assumed  to  result  from  the  viscosity  of  the  metallic  jet 
due  to  the  velocity  gradient  within  the  jet  and  the  resulting  shear 
stress  within  the  jet.  Also,  the  Reynolds  numbers  defined  by  the  two 
models  are  identical  for  the  penetration  process  when  the  viscosity  of 
the  target  is  the  pertinent  parameter. 

The  radial  momentum  equation  is: 


3t 


3r 


3z 


1 fl  3 , 

— Z ;=  Cr  o ) 

u r 3r 

o L 


and  for  Re  = 0(1/R  ) and  k = 0(R),  with 


•^rr  . 2 3V 

1 ^ k Re 

p u ^ r r Re  3r 

0 

J.  ^ 

/2 

&'  ■ ’ 

f'sv.^y  VsuV 

(3z  3r/  \3z/ 

1/2J 

20 


3V 

Xi  3V 

- 3V 

— + 

V — 

+ u — 

at" 

3r 

3z 

using  the  continuity  equation. 
For  k « O(R^) . 


at  3r  3z  ® j_r  3r  3r^  J ® |_3r^  r 3zJ 


Thus,  the  continuity  and  momentum  equations  for  a visco-plastic  model 
reduce  to  the  equations  for  a Newtonian  fluid_for  small  values  of  the 
normalized  yield  stress  in  pure  shear,  i.e.,  k = O(R^),  where  the 
Reynolds  number  may  differ  between  the  two  models,  as  discussed  earlier. 


VII.  DISCUSSION  OF  EQUATIONS 

The  governing  equations  formulated  for  the  Newtonian  fluid  model 
and  for  the  visco-plastic  model  have  been  simplified  by  assuming 
incompressible  flow,  axisyrametric  flow,  and  by  retaining  only  dominant 
order  terms.  However,  the  final  equations  are  still  nonlinear  and 
analytical  solutions,  even  approximate  solutions,  are  not  readily 
apparent . 

Nevertheless,  the  governing  equations  for  the  Newtonian  fluid  are 
equivalent  to  the  transient  boundary  layer  equations.  The  boundary 
layer  equations  have  been  treated  extensively  in  the  literature.  How- 
ever, the  boundary  and  initial  conditions  differ,  of  course,  between 
shaped  charge  jets  and  boundary  layer  flow.  For  steady  state  conditions, 
the  equations  can  be  solved  by  similarity  transformations  for  generalized 
boundary  conditions.  Steady  state  solutions  however,  may  not  realisti- 
cally describe  the  true  jet  behavior. 

The  visco-pla.stic  model  governing  equations  are  similar  to  the 

Newtonian  fluid  model  except  that  an  additional  term  k (for  IT  * 0(R))is 
present.  r 

If  a one  dimensional,  transient  flow  is  assumed  and  the  normal  stress 
terms  are  retained,  a somewhat  simpler  set  of  equations  results. 

Chou  and  Carleone**’ are  investigating  the  solution  of  the  one  dimensional 
transient,  inviscid,  incompressible  equations  using  a strain  hardening 
model  for  the  normal  (axial)  stress.  Their  equations  employ  a 
Lagrangian  coordinate  system. 

^^Chout  P.  C,  and  J,  Carteone,  "The  Stability  and  Breakup  of  Shaped 
Charge  Jete,"  Submitted  to  Third  International  Sympoeium  on  Ballietiaej 
Karlarue,  Germany,  March  1977. 
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VIII.  CONCLUSIONS 


Governing  equations  have  been  developed  for  an  Eulerian,  axi- 
synmetric,  incompressible  flow.  The  stress  tensor  was  modeled  both  as 
a Newtonian  fluid  and  as  a visco-plastic  material. 

The  visco-plastic  material  governing  equations  were  shown  to  relax 
to  the  Newtonian  fluid  equations  when  the  yi'ild  stress  in  shear  was 
assumed  very  small.  However,  the  Reynolds  numbeis  defined  in  the  two 
models  are  not  necessarily  identical. 

For  both  models  appropriate  nondimensional  groups  were  defined  and 
the  dominant  terms  were  determined  by  an  order  of  magnitude  analysis. 
Other  models  of  the  stress  tensor  can  be  analyzed  by  a similar  order  of 
magnitude  analysis  as  long  as  the  stress  tensor  can  be  expressed  as 
functions  of  the  velocities  or  velocity  gradients,  (or  displacements). 

The  assumption  of  incompressible  flow  is  extremely  important.  This 
assumption  decouples  the  continuity  and  momentum  equations  from  the 
energy  and  state  equations.  The  energy  equation  given  for  the  Newtonian 
fluid,  incompressible  flow  model  involves  an  internal  heat  generation 
term.  Terms  of  this  nature  (frictional  heat)  are  not  included  in  the 
energy  equations  formulated  for  the  hydrocodes.  This  term  will  include 
the  jet  temperature  and  may  influence  the  balance  between  the  kinetic 
energy  and  the  internal  energy. 
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APPENDIX  A 


THE  JET  REYNOLDS  NUMBER 

The  derivation  of  the  final  momentum  and  energy  equations  requires 
that  the  Reynolds  number  be  relatively  large  to  allow  the  inclusion  of 
the  viscous  terms  in  the  final  order-of-magnitude  equations.  This 
assumption  is  now  investigated. 

The  Reynolds  number  is  given  by 

p u L 

Re  = — : 

u 

a representative  jet  tip  velocity  may  range  from  3 - Dnon/usec, 

a representative  jet  length  (measured  from  the  virtual  orgin)  may 
range  from  100  to  600mm, 

a representative  kinematic  viscosity,  p/p  can  range  from  0.25  to 
5,58,9»i**  inin^/ysec  for  most  jet  or  target  metals  of  interest. 

Thus,  the  Reynolds  number  may  range  from  10^  to  10^.  Usually  a Reynolds 
number  greater  than  10^  should  suffice  for  a large  Reynolds  number 
assumption. 


Harlow  and  W.  Praaht,  "Formation  and  Penetration  of  High-Speed 
Coll-apee  Jets,"  Fhgeiae  of  Fluids,  Vol,  9,  Ho,  10,  October  1966, 
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